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INTRODUCTION 


Recently,  much  has  been  done  on  the  flow  problems  where 
chemical  reactions  take  place  in  the  gaseous  mixture. In 
these  types  of  flows,  one  must  consider  the  chemical  reactions 
as  well  as  the  fluid  dynamics  simultaneously.  Clarke^  has  ob¬ 
tained  the  linearized  solution  for  a  flow  past  a  corner  in  a 
reacting  gas  mixture ,  while  Vincent!^  solved  the  problem  of 
nonequilibrium  flow  over  a  wavy  wall  and  showed  that  the  extra 
drag  due  to  nonequilibrium  effects  is  positive.  Li  and  Wang^ 
have  presented  a  linearized  slender-body  theory  of  dissociating 
gas  flows.  But  no  numerical  results  had  been  worked  out.  This 
paper  will  be  concerned  with  some  quantitative  results  of  linear 
slender-body  theory.  In  particular,  the  drag  of  two  cases  will 
be  considered;  namely,  a  circular  cone  and  a  parabolic  spindle. 
These  results  should  be  useful  for  the  development  of  theories 
with  higher  approximations.  In  fact  the  extra  drag  due  to  non¬ 
equilibrium  effects  turns  out  to  be  quite  appreciable  in  the 
first  case.  For  the  case  of  the  circular  cone,  the  extra  drag 
is  positive,  while  for  the  case  of  the  parabolic  spindle,  the 
extra  drag  is  negative  indicating  that  the  nonequilibrium  effect 
decrease  the  total  drag. 


GOVERNING  EQUATIONS  AND  BOUNDARY  CONDITIONS 


In  the  following  we  shall  be  concerned  with  flows  where 
the  only  chemical  reactions  taking  place  are  molecular  dis¬ 
sociations  and  atomic  recombinations.,  By  ruling  out  ionization, 
we  are  actually  restricting  ourselves  to  a  certain  temperature 
range. (10)  Furthermore,  we  assume  that  viscosity,  diffusion  and 
heat  conduction  of  the  gases  are  negligible.  From  this,  the  con 
tinuity  equation  for  the  mixture  is 


where  P  is  the  density  of  the  mixture  and  u.  is  the  macroscopic 
velocity  vector  of  the  mixture. 


The  dynamical  equation  of  motion  is 


= 


(2) 


Here  is  the  substantial  derivative  and  P  is  the  pressure. 


The  energy  equation  for  the  mixture  becomes 


P  Dh  DP_ 

'  Dt  Dt  (3) 


where  f)  is  the  specific  enthalpy  of  the  mixture.  From  the  dyna 
mical  equation  of  motion  and  the  expression 


H  =  ±  U.-0L  +  h 


(4) 


where  H  is  defined  as  the  stagnation  enthalpy,  the  above  energy 
equation  can  also  be  written  as 
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P^H_=  JP 
'  Dt 


(5) 


From  the  thermodynamic  point  of  view,  we  asstime  the  different 
gases  in  the  mixture  to  be  thermally  perfect.  Thus  we  can  write 
the  thermal  equation  of  state  for  a  dissociating  gas  as 

P  =  (i+c*')  f  T  (6) 

O 

where  is  the  molecular  gas  constant,  T  is  the  temperature 
of  the  mixture  and  o<  is  the  degree  of  dissociation  defined  as 
the  mass  fraction  of  gas  dissociated.  The  caloric  equation  of 
state  is 


e  =  e  o  jT) 


(7) 


where  6  ,  the  specific  internal  energy  of  the  mixture,  consists 
of  the  specific  internal  energies  of  the  participating  species. 
For  dissociating  diatomic  gases, 

e  =  (i-*)e2  +  (8> 

where  6^  and  0^  are  the  specific  internal  energy  of  the  mole¬ 
cular  and  atomic  species  respectively.  With  the  assumption  that 
the  translational,  rotational  and  vibrational  degrees  of  freedom 
are  fully  excited,  and  the  contributions  of  energy  from  the  other 


4 


degrees  of  freedom  are  negligible,  we  can  express  the  specific 
internal  energies  derived  from  statistical  thermodynamical  con¬ 
siderations  as. 


£it  =  +  +\'CX^ 


eu-  R; T 


=•  0e_|  RiT  “  f(9)R:T 


(9) 

(10) 
(ID 


where 


d  •  dissociation  energy  per  unit  mass  of  diatomic  gas 
R  »  universal  gas  constant 
R.  ss  gas  constant  per  unit  mass  of  i-species 
W; 551  molecular  weight  of  i-species 


e  = 


hV 

fcT 

e 


f(e)=  -^rrf 

^  »•  Planck's  constant 

k*  Boltzman's  gas  constant 

^  =  vibrational  frequency 

T  “  1  when  i  =  | 

1  *=  0  when  i  = 

and  the  subscript  i  may  stand  for  /  or  for 
dissociating  diatomic  gas.  In  this  case,  the 
of  state  is 


the  case  of  a 
caloric  equation 


e  =  [-f-+  fieui-«>+A-]-^T  +  o<G( 


(12) 
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The  specific  enthalpy  is  expressed  by 

A  =[ir +  (13) 

In  the  case  of  ideal  dissociating  gas,  f  (*)»■£ 

The  continuity  equation  for  the  i-species  can  be  written 
as 

~^jr  +  V  •  (  ?;  ui)  =  CT-  (1^) 

where  0“£  is  the  rate  of  mass  product  of  i-species  per  unit 
volume,  or  usually  known  as  the  chemical  source  term  of  i-species 
and  IT-  is  the  macroscopic  velocity  of  the  i-species.  For  ex¬ 
ample,  in  a  dissociating  gas,  the  conservation  of  atomic  gas  can 
be  expressed  as 

+  V-[fcX  (UfVA)J=  A  (15) 


where  /)  is  the  atomic  mass  rate  production  per  unit  volume,  and 
V*  is  the  diffusion  velocity  vector  of  the  atoms.  In  regions 
outside  the  boundary  layer  or  where  the  flow  is  inviscid,  the 
term  involving  ^  can  be  neglected.  Hence 


2*  .  JL 
Dt  ^ 


(15a) 


From  consideration  of  chemical  kinetics,  the  term  on  the  right 
hand  side  of  the  equation  can  be  obtained.  Following  J.F.  Clarke^, 
we  write,  for  the  case  of  a  diatomic,  dissociating  gas 
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P2L-4KvfVgQ  vk  ,  ,  Cl  i  f  .  ,1 

vt  Wj.  [k^  y[k0'^“-  * j 


(16) 


where  Kj  and  K,*,  are  dissociation  and  recombination  rate 


constants  respectively,,  Also, 


T  = 


WC 


fl(i+oO 


(17) 


known  as  the  characteristic  reaction  time,  and 


K  = 


Kf  iVt 

4f 


(18) 


which  can  be  related  to  the  atom  mass  fraction  under  equilibrium 
conditions  in  some  circumstances. 


If  we  restrict  ourselves  to  steady  flow,  the  above  governing 
equations  for  inviscid,  dissociating,  diatomic  gas  with  the  stated 
assumptions  become 


V*(f  Z)  «  o 

(19) 

(ZL- v)TZ  -  -  J  vp 

(20) 

7L-V  H  =  0 

(21) 

P=o+«)f-^r 

(22) 

A  =  h  (Tj  cx) 

(23) 

-  .  V(y  =  _L  j\o-<x)-<XlJ 

(24) 
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The  specific  entropy,  s,  can  be  introduced  by  the  equation 


7”  ds  =  d  /i  “  |3  }  do( 


(25) 


where  is  the  difference  of  the  chemical  potentials  of 

the  molecules  and  atoms. 


'yU-z,  can  be  related  to  o<  and  T  as  4,5 


These  equations  can  be  combined  to  yield^ 


-4  V-7t  -  •  ycx 

where 


is  the  frozen  speed  of  sound,  and 


In  the  above  expression 


(27) 


(28) 


(29) 


(30) 


is  the  specific  heat  at  constant  pressure  and  frozen  composition 


and 


s  =  -  J-  (2L) 

n  tvvrM 


(31) 
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is  the  volume  expansion  coefficient  at  constant  pressure  and 
frozen  composition.  The  problem  now  is  to  solve  equations 
(24)  and  (27)  simultaneously,  subject  to  the  proper  boundary 
conditions;  namely,  that  the  flow  be  tangent  to  the  body  on 
the  surface  of  the  body  and  that  there  be  no  chemical  reaction 
taking  place  between  the  gas  and  the  body. 


Linearization 

For  flow  past  slender  bodies,  the  disturbances  are  very 
small.  To  simplify  the  problem  of  steady  flow  of  a  simple 
dissociating  gas  mathematically,  we  may  assume  the  disturbances 
to  a  free  stream  in  chemical  equilibrium  to  be  small  and  pro¬ 
ceed  to  linearize  the  equations.  We  assume  that  not  only  will 
the  perturbation  velocities  be  small,  but  the  deviations  of 

«  ,  from  their  respective  free  stream  values  will 

/  /  * 

also  be  small.  Let  u  ,  \J  and  be  the  non-dimensional  per¬ 
turbation  velocity  potential  defined  as 
'  u. 


UL  - 


U. 


i r  = 


i r 

LAoo 

uS 


(32) 


In  the  linearization,  we  neglect  terms  involving  squares  or 
products  of  the  perturbation  quantities.  The  validity  of  these 
assumptions  holds  for  slender  bodies.  With  these  approximations. 
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equations  (24)  and  (27)  can  be  combined  to  give*-! 

+  IK#  +  +  « 


(33) 


where 


Hoc 


Me*  ~ 

?  - 
7  - 
5  - 


and 


(34) 


(35) 


stands  for  the  ratio  of  characteristic  flow  time  and  chemical 
relaxation  time.  (_  is  a  suitable  characteristic  length  of  the 
flow.  Within  the  linearized  theory, ^  equation  (33)  can  be 
written  as 


(36) 


where  the  non-dimensional  velocity  potential  &  has  been  introduced. 
The  non-dimensional  perturbation  velocities  are  given  as 


(37) 
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The  solution  of  equation  (36)  subject  to  the  linearized  boundary 
conditions  has  been  solved  by  Li  and  Wang^  by  using  the  method  of 
Laplace  and  Fourier  Transforms. 

Drag 

From  the  expression  for  , the  perturbation  velocities, 
surface  pressure  and  hence  the  drag  coefficient  for  slender  bodies 
can  be  obtained  as^ 


Where  S(x)  denotes  the  cross  sectional  area  of  the  body  and  should 
cause  no  confusion  with  the  specific  entropy.  The  first  four  terms 
are  the  same  as  those  obtained  by  Ward  and  the  last  two  terms  are 
purely  due  to  non-equilibrium  effect.  The  third  term  can  be  ob¬ 
tained  when  specific  body  shapes  are  considered.  In  particular,  for 
bodies  of  revolution  and  omitting  the  base  drag  for  bodies  not 
pointed  at  both  ends,  the  expression  for  the  drag^  becomes 
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if.x  ^ 


f  jes'(x)  5  "(t)  >*- 5 1  jf  Jx- 


•  0  v/r*  ,  p-*<*-?>_  .>-#.<*-}> 

+  «4Srf)  +  if/  /  S'^'W- - — - ^ 


o  e 


:'t 


.swjS'd)- - — - ^ 


(39) 


where  o<0  is  the  angle  of  attack,  and 

a  - 

r1  £.  U„  /w£ 

Mu  C  *£ 

P  =  Xc  tf  „  ML  L  Mfi  -  /  _ 

The  last  two  terras  in  the  above  expression  (Eq.  39)  gives  the 
extra  drag  due  to  nonequilibrium  effects.  The  relative  importance 
of  the  drag  will  be  investigated  for  two  cases;  namely,  the  unyawed 
circular  cone  and  the  parabolic  spindle. 


Examples 


In  the  following  examples,  computations  are  based  on  the 


va  lue  s 


_z  E>e 

£>  =  o.  *  =  1.5  and 


f  ^fbo 

l  % 


=  1.35.  Correspond 


to  these  values,  the  equilibrium  and  frozen  Mach  numbers  are  2.12 


and  1.83  respectively. 


Case  a)  Unyawed  cone 
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The  profile  of  the  cone  is  given  by  (see  figure  1) 
t<%ri  &  =  (40) 

where  &  is  the  semi-vertex  angle  of  the  cone.  Hence, 
zS  (*)  =  TTR2--  TTx" 

Knowing  S '  (x)  and  S"(x)  the  drag  of  a  cone  in  a  simple  dis¬ 
sociating  gas  can  be  calculated  from  equation  (39)  by  double 
integration.  For  this  case,  denoting  and  d(2)  as  the  drag 

of  the  cone  in  an  inert  gas  and  the  extra  drag  due  to  chemical 

Da> 

reaction,  the  expression  for  TX  IT*"-  turns  out  to  be 


X 

X 


D  T"  ~  **  j 

^  -6  . 

1  i  4  i~e  ( 1 )  _  i-  e  "(fil-H) 


+[(P/+7fb}]^'"e"P  w^'"e  P 


(42) 


The  expression  for 

z 

result  obtained  as 


can  be  also  integrated  and  the 


— - — Jr  =  &  '  i  “  £"  ~ ~) 

i  f-  L  V  - 


(43) 


M U 


Too  ^mo 


is  essentially  the  inverse  characteristic 


PL- 

relaxation  distance.  For  a  model  with  length  A  equal  to  the 
characteristic  relaxation  distance,  i.e.,  >  the  ratio 
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has  been  calculated  for  *  5°  to  30°  at  5°  inter¬ 


vals  c  The  result  is  given  in  table  10  For  a  model  of  fixed 


length,  the  effect  of  chemical  reaction  is  seen  to  decrease 
rapidly  for  increasing  .  (Table  2). 


Case  b)  Parabolic  spindle  at  zero  angle  of  attack  (Fig.  2). 
The  profile  of  this  body  is  given  by  the  equation 


R<*)  =  j) 

Hence 

SW-  7TR  -  7 riV(/-y)z 


(44) 


(45) 


From  the  expression  for  drag  given  by  equation  (39),  and  noting 
that  S'(l)  =  0, 


(46) 
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and 


(47) 


(« /  A 

The  ratio  of  D / $'  for  A. a*  is  shown  in  table  3.  In 

this  case  this  ratio  is  independent  of  the  slenderness  ratio. 


Conclusion 

(2) 

For  the  case  of  the  cone,  Dv  '  turns  out  to  be  positive 

and  is  quite  appreciable  compared  to  d(^ .  This  seems  to  be 

12 

in  support  with  the  reasoning  in  Kusukawa  and  Li's  report; 
namely,  the  presence  of  chemical  reaction  will  increase  the  en- 
tropy  just  like  in  the  case  of  vorticity.  The  ratio  /d 
for  BJl  =  1  increases  quite  rapidly  for  increasing  semi-vertex 
angle  of  the  cone.  In  particular,  the  value  already  reaches  107, 
for  a  cone  of  semi-vertex  angle  of  between  15°  and  20° „  This 
result  justifies  the  investigation  of  theories  of  non-equilibrium 
flow  with  higher  approximations.  These  higher  order  theories 
are  useful  not  only  in  the  increase  of  accuracy,  but  also  in 
their  applicability  to  a  wider  range  of  problems.  For  example, 
by  the  higher  order  theories,  we  are  able  to  extend  the  Mach 
number  range  so  as  to  approach  the  hypersonic  case.  The  effects 
of  chemical  reactions  decreases  very  rapidly  for  a  given  model 
if  |3  is  increased  as  shown  in  table  2. 


15 


i 

l 

I 

1 


1 


In  the  case  of  a  parabolic  spindle,  d(^)  is  negative  and 
is  small  compared  to  From  thermodynamic  considerations, 

it  has  been  shown  in  Ref.  12  that  the  extra  drag  caused  by  chemi¬ 
cal  equilibrium  will  be  positive  for  a  subsonic  flow  past  a  body. 
In  order  to  relate  the  subsonic  and  the  supersonic  cases  for 
slender  pointed  bodies,  the  following  expression  for  the  extra 


drag  obtained  from  reference  9  can  be  used 


tr* 


3? I  ISMS® 


c  0 


dr  dx_ 


(48) 


for  both  supersonic  and  subsonic  case.  Since  Pi<  fit  for  sub¬ 
sonic  case  and  for  supersonic  case,  the  sign  within  the 

bracket  is  reversed  for  the  two  cases.  Hence,  for  the  same  body, 
if  the  subsonic  drag  due  to  chemical  reactions  is  positive,  its 
counterpart  in  the  supersonic  case  must  be  negative.  It  was 
shown'*''*'  that  in  a  medium  where  chemical  reaction  is  taking  place, 
disturbances  are  being  damped.  Hence,  in  a  supersonic  flow,  the 
wave  drag  in  flow  of  a  reacting  gas  past  a  body  will  be  less  than 
that  of  a  non-reacting  gas.  As  pointed  out  at  the  end  of  refer¬ 
ence  12,  this  phenomena  would  have  to  be  included  in  the  dis¬ 
cussion  of  the  extra  drag  from  effects  of  chemical  reaction.  The 
results  on  the  parabolic  spindle  indicates  that  the  decrease  of 
wave  drag  is  sufficiently  large  to  cause  D^)  to  be  negative.  Un¬ 
like  in  the  case  of  the  cone,  where  only  dissociation  takes  place, 
in  the  case  of  parabolic  spindle  both  dissociation  and  recombination 
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occurs.  In  addition  to  this  difference,  the  ratio  D  /d° 
in  the  latter  case  is  independent  of  the  slenderness  ratio.  For 
a  fixed  ,/  ,  as  in  the  case  of  the  cone,  the  effects  due  to 

chemical  reaction  decreases  rapidly  for  increasing  ^  ,  which 
may  be  interpreted  as  the  inverse  characteristic  relaxation  dis¬ 
tance.  The  calculations  in  this  paper  are  based  on  1.83 

and  M.*,  =  2.12  so  that  the  linearized  theory  is  applicable.  The 
results  as  indicated  here  may  be  quite  different  from  the  hypersonic 
slender  body  behavior.  For  that  purpose,  higher  order  theories 
with  a  higher  Mach  number  range  of  applicability  must  await  future 
development. 


>1 


Parabolic  Spindle 
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Table  1:  Variation  of  d(2)/d(1)  with  &  for  an  unyawed  cone. 

I  /=  IS  o* 


<5 

5° 

h* 

o 

o 

15° 

20° 

25° 

30° 

4.3 

6.2 

8.6 

12.1 

17.7 

29.7 

Table  2:  Variation  of  D^2VD(1)  with  ^  for  an  unyawed  cone, 

15  S  =  5° 


ft 

1 

15 

.0715 

.715 

7.15 

DU>  .  7 
tf7"1  / 

4.3 

4.0 

.11 

.11  x  10“2 

Ratio  of  D^2Vd(1) 
for  a 

Table  3  parabolic  spindle, 

A/-i 


(Z) 

D 


_<o 


.000109 


=  .14% 


.081 
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SYMBOLS  USED 

0.^  frozen  speed  of  sound 
Qe  equilibrium  speed  of  sound 
f\  atomic  mass  rate  production  per  unit  volume 
Cp^  specific  heat  at  constant  pressure  and  frozen  composition 
D  •  drag 

D'j  drag  of  flow  past  a  body  in  inert  gas 
£)J)  extra  drag  due  to  chemical  reaction 

4  dissociation  energy  per  unit  mass  of  diatomic  gas 
Q  specific  internal  energy  of  the  mixture 

5  specific  internal  energy  of  the  molecular  species 
2. 

e/  specific  internal  energy  of  the  atomic  species 
h  specific  enthalpy  of  the  mixture 
k)  Planck's  constant 
H  stagnation  enthalpy 
k  Boltzman's  constant 
Kf  dissociation  rate  constant 
K.  recombination  rate  constant 
L  characteristic  length  in  the  flow 
M,  frozen  free-stream  Mach  number 

Too 

equilibrium  free-stream  Mach  number 
P  pressure  of  the  mixture 
P  universal  gas  constant 

p“.  gas  constant  per  unit  mass  of  i  -species 
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Z >  when  i  ~  z 


A 

chemical  potential  of 

atoms 

A* 

chemical 

potential  of 

molecules 

2. 

B  = 

Bte1 

■ 

£>e. — 

R7" 

i”  .  Bt=  J 

3  = 

< 

r  1 1 

A  _ 

V 

L  m£-  1  J 

1  ' 

s  - - L2 - 

Me*! 

f(*)= 

i>  _  J5i  Wt 

N  ~  Kv'4f 


5  =  ^  ,  ;  -  #/l  ,  5  =  yL 

SUBSCRIPTS 


t  translational 
/V  rotational 
^  vibrational 
{  frozen 
€.  equilibrium 
free-stream 
Z  molecular  species 
I  atomic  species 
i  i-species 
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